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Highly accurate methods such as coupled cluster (CC) techniques can be used for periodic systems 
within the framework of the method of increments. Its extension to low-dimensional conducting 
system is considered. To demonstrate the presented approach a clean Mg (0001) surface is selected, 
where the CC treatment with single and double excitations and perturbative triples is used for 
calculation of the surface energy. Further example concerns the adsorption energy of Xe on the 
metal surface. The obtained results can be used to verify performance of the approximate methods. 
Along with the computational speed-up at the high-level of accuracy, application of the method of 
increments provides for a possibility to analyze influence of individual correlation-energy increments 
on the studied property. 



Although density functional theory (DFT) remains the most frequently used method in the field of solid-state 
physics and material science, there is no clear way to systematically improve results obtained with the present-day 
density functionals. This yields preferences to the wavefunction-based correlation techniques, a hierarchy of methods 
that, starting from the one-electron Hartree-Fock (HF) solution, allows one to approximate the many-body nature of 
the true ground-state wavefunction in a rigorous and systematic manner. One of the largest problems standing on the 
way of the application of such methods to extended systems is their unfavorable scaling with respect to the increasing 
system size. From a physical point of view, however, this difficulty should be avoidable because the correlation hole 
around an electron is a fairly local object. 

The reformulation of the many-body wavefunction in terms of localized orbitals defines the group of the so-called 
local correlation methods (see e.g. Refs. [1-7 ). One method of this type, the method of increments [7 , combines 
HF calculations for periodic systems (which give Coulomb and exchange terms of the total energy) with correlation 
calculations for finite embedded fragments of a solid, and the total correlation energy is written as a cumulant 
expansion in units of contributions from localized orbital groups of increasing size: 



The summation over i includes all orbitals in the reference cell, while j and k extend over the entire system. The 
one-body increment (e^) involves all excitations from the orbitals while the two-body increment (Ae^j) includes 
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the nonadditive excitations summed over orbitals i and j. In order to get the latter quantity, one subtracts from 
the correlation energy, obtained for the joint orbital system zj, the corresponding one-body contributions: Asij = 
Sij — {si -\- Sj)' Higher-order increments are defined in analogous way, e.g. for the three-body term we get: Asijk = 
^ijk — i^i +^/e) — i^^ij + ^^ik + ^^jk)- the limit of the sum over all z, j, k as well as over orders of increments, 
one would obtain the exact correlation energy per unit cell of the infinite system. In practice, the neglect of four-body 
and higher-order contributions is possible and the incremental expansion converges swiftly with distance (see Ref. [8] 
and the references therein). The translation symmetry of crystals reduces the amount of increments to be evaluated 
drastically. In order to use this advantage, one determines the so-called weight factors {uj) with which the individual 
increments occur in the considered structure. 

The presence of well-localized orbitals is a precondition for the applicability of local correlation schemes. For metals, 
it is no longer possible to generate well-localized Wannier orbitals by unitary transformation within the occupied HF 
space. To overcome this difficulty it has been suggested to base the incremental expansion on a well-localizable model 
system rather than on the real metal and to allow for delocalization (and thereby approaching the real metal) only 
gradually within the various levels of the many-body expansion [9] . This scheme is tested for series of bulk metals 
and the obtained cohesive properties are shown to be in very good agreement with experimental results p!QHT4] . 

Due to the need to calculate the correlation energy increments in finite fragments of the solid, one has to construct 
an embedding, which simulates the conditions in the infinite system. The embedding is chosen in such a way that the 
localized orbitals in the finite systems are quite similar to the Wannier orbitals of the periodic system. Type of the 
embedding depends on the nature of binding in the system under consideration, e.g. dangling bonds saturation with 
hydrogen atoms in the case of covalent systems (see e.g. Refs. [l5l[T6]) or a point charge array in the case of an ionic 
crystal (see e.g. Refs. [T71 [13 )• When dealing with metal, the solution is not obvious since the correlation effects 
that inffuence the ground-state properties are quite local in this case. However, an attempt to get good agreement 
with experiment for Mg-bulk correlation energy using the free clusters fails due to the mixing of species of different 
nature [11 . That is way different variants of the embedding and the localized orbitals used for the incremental scheme 
have been developed and tested; their advantages and disadvantages are discussed in Refs. [9( [14]. 

Being defined for three-dimensional case, the method of increments can nevertheless be applied for systems with 
reduced symmetry. For instance, to quantify the molecule-surface interaction the adsorption energy can be written as 
£^ads = Es-\-A — Es — Ea = ^id^ + ^ads^' whcrc Es-\-A is the total energy of the substrate/adsorbate system, and Es 
and Ea are the energies of the fragments at the same coordinates as in the substrate/adsorbate system corrected for 
the basis set superposition error. Again, the HF energy, is calculated for the periodic system in the standard 

way, and the correlation contribution, E^^f ^ is calculated within the incremental expansion: 



Here: r]A is the correlation contribution to the adsorption energy within the adsorbate, r]i and r]ij are the changes in 
the surface increments due to adsorption, rjAi are caused by the interaction between localized orbitals of the adsorbate 
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and the surface. Overall, one has to take into account all orbital groups, that change due to the interaction of an 
adsorbate with a substrate. This approach can be especially efficient in the cases where DFT is less successful, e.g. 
when dealing with physisorption [161 EHl IE] , although there is no doubt that chemisorption can be investigated as 
well. 

Combining the incremental approach with the scheme proposed for metals, it should be possible to consider metallic 
surfaces. In this manuscript, we present the results on the Mg(OOOl) surface and the Xe/Mg(0001) system, where the 
above approach is used to investigate surface stability as well as interaction energy between the metal surface and the 
adatom. 

II. COMPUTATIONAL DETAILS 

The periodic mean-field calculations for the studied system were performed with the program package CRYS- 
TAL 09 [28 . In order to obtain converged results for the HF binding energy, the default parameters were modified. 
For the pseudo-overlap tolerances IT0L4 and IT0L5, used for the prescreening of the exchange integrals, tighter 
values than usual have been employed: 15 and 80, respectively. The other CRYSTAL tolerances (ITOLl-3), used 
for the screening of Coulomb integrals, have been set to 7. The shrinking factors defining the /c-mesh were set to 16 
and 32. The basis sets employed are the crystal-optimized all-electron basis set of valence-double-^ (VDZ) quality 
including polarization functions [29] for Mg and the one optimized for the 8- valence electron pseudopotential [30 for 
the Xe atom. The adsorbate concentration is set to approximately 3.7-10"^ A~^. The converged result for the surface 
energy at the HF level is achieved for the slab consisting of feefteen Mg-layers. When calculating adsorption energy, 
the seven-layer Mg-slab was used. 

For compatibility, the DFT calculations reported here were performed with the same thresholds and basis set as used 
for the HF calculations. For the exchange-correlation energy functional we employed the local density approximation 
(LDA) [24| and the generalized gradient approximation as parameterized by Perdew et al. (PBE) [25 . In the context 
of DFT, the contributions of the long-range van der Waals interactions on the adsorption energy of the xenon on 
magnesium were accounted for by means of the semi-empirical approach proposed by Grimme [27 . 

The correlation-energy increments are calculated for selected fragments which reflect the geometry of the crystal. 
These fragments have two components, firstly the atoms i, j, /c, ... to be correlated (in the center of the fragment) 
and secondly the embedding atoms. We select the positions of the atoms to be correlated and surround each of 
these atoms by embedding shell with a cutoff for the embedding of 2 times the in-plane lattice parameter of the 
hep Mg. This results in the cluster consisting of 51 atoms (1 atom in the central part and 50 embedding atoms), 
when considering the one-body contribution. To get the results presented in this work the embedding was partitioned 
into two regions, a true embedding region (second embedding shell), where the atoms are described with a minimal 
basis, and an intermediate region (first embedding shell), where the occupied orbitals are kept frozen, but additional 
basis functions are supplied for the re-optimization of the orbitals in the central part of cluster. An example of the 
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embedded two-body cluster is shown in Fig. [l] 

Correlation calculations are performed with the program package MOLPRO [ST, using the coupled-cluster treat- 
ment with single and double excitations and perturbative triples [CCSD(T)] [32]. For the perturbative approaches in 
the incremental scheme, we have re-canonicized the space of the localized orbitals which are correlated to obtain a 
diagonal Fock matrix. The basis sets employed are the cc-pVTZ of Dunning [33] for Mg and the cc-pVTZ optimized 
for the 8-valence electron pseudopotential [34 increased by even-tempered basis functions to [454p3(i2/] for the Xe 
atom. Embedding atoms are described with the 2-valence electron pseudopotential [35 and either the corresponding 
basis set of VDZ quality (first embedding shell, see Fig.jl]) or a minimal basis set (second embedding shell, see Fig.jl]). 

III. RESULTS AND DISCUSSION 

Surface energy (£^surf) quantifies the disruption of bonds that occurs when a surface is created. Within the slab 
approach, mathematically this can be written as follows: £^surf = (£^siab — ^ x £^buik)/2 (n is a number of layers in the 
slab). Therefore in the formalism of the method of increments, correlation contribution to the surface energy (7) can 
be determined as a difference between the increment calculated using the "surface" cluster (Fig. [T]), multiplied with 
its weight factor (cj®^^^), and the corresponding one calculated using the "bulk" cluster, multiplied with its weight 
factor (cj^^^^). For example, for the two-body correlation-energy increment one writes: 

A surf surf a bulk bulk /o\ 

-iij Asij X ujij - Asij X ujij . (3) 

In the physics of solids, surfaces must be intrinsically less energetically favorable than the bulk of a material. From 
the computational point of view this means positive sign for the value of the surface energy. The experimental result 
for the Mg(OOOl) surface is +0.383 eV/atom [20]. HF method, neglecting electron correlation effects, yields only 47% 
of Esuif (Tab.[l|. Summing up 1-, 2-, and 3-body contributions one gets E^^l"^ = +0.20361 eV (Fig. [2|, that leads to 
^sSJf = +0-385 eV (Tab.|l|, that is in very good agreement with the experimental value. 

The adsorption of rare-gas atoms on metal surfaces represents prototypical physical adsorption systems and, thus, 
Xe adatom on Mg(OOOl) can be viewed as a good model to test our approach. Two adsorption sites have been 
investigated: (i) a low-coordinated atop position and (ii) the fee hollow site (The corresponding dissociation curves 
are presented in Supplementary Fig. SI). The obtained results are summarized in Tab. [ij Whereas HF shows 
no binding in both cases, correlations stabilize the structures, giving preference to the atop geometry. Up to our 
knowledge, there are no experimental data regarding the studied Xe/Mg(0001) system, however, the atop position 
was experimentally found to be more favorable, when investigating Xe adatom on Cu(lll), Pd(lll), Pt(lll), and 
Nb(llO) [2T1122]. 

It is interesting to compare the results evaluated by means of the presented approach with those of DFT. The 
data obtained employing LDA and PBE are listed in Tab. [l| As usual, LDA overestimates the binding, that leads 
to too short lattice parameters of the hep-Mg and too large ^coh- Further, surface is c.a. 4% less stable than it 
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is experimentally and £^ads is strongly overestimated compared to the benchmark. Although PBE results are in 
reasonable agreement with experiment for bulk characteristics, this functional yields overestimation of the surface 
stability, demonstrating behavior typical for gradient-corrected functionals [26]. Finally, PBE fails to describe Xe 
physisorption, that is not surprising as well. In the latter case, performance of PBE can be improved when employing 
DFT-D2 scheme of Grimme [27] . This way, adsorption energies and equilibrium distances are in good agreement with 
the presented CCSD(T) results, deviating by less than 1 %. 

An asset of the incremental scheme is that along with the computational speed-up at the high-level of accuracy one 
gets the possibility to analyze various contributions to the correlation energy. This gives us a chance to penetrate 
deeper into the nature of binding and answer some interesting questions. The first one is: why do correlations 
destabilize the surface? One can see from Fig. |2j that E^^^^ is mainly determined by the two-body contributions. 
Therefore, in order to find an answer to this question, it is reasonable to analyze the two-body correlation-energy 
increments considered in this work (Tab.[ll|). Indeed, Ae^^^ are higher in magnitude as compared to the corresponding 
contributions of bulk. At the same time, in most cases is larger than 0, that is due to the lower weight factors 
ct;|^^^ Weight factor for the two-body increment is defined as the number of neighboring atoms located at the certain 
distance from the one set as a central (divided by two in order to get the contribution per one atom). In Fig. [l] 
where the cluster used for the calculation of Aei2 is presented, the six distances of one of the central Mg in the 
"bulk" cluster are shown with thicker lines; three of them, drawn with dashed lines, will disappear, when constructing 
the corresponding "surface" cluster. As a result, the corresponding weight factor is reduced by 2. The situation is 
different, when considering "surface" clusters with two atoms lying in the top, (S), layer. In such cases (footnoted 
with a in Tab. [n|), u^^^^ = ujf-^^^ leading to negative jij. Summing up all 7, one gets E^^l^^ which is large than 0. 

The second question can be formulated as follows: Noting that Xe adatom has closed shell, it can be assumed that 
this adsorbate should prefer maximally coordinated sites (minimizing the Pauli repulsion). What is the reason for the 
discrepancy between this expectation and our observations? In order to understand this phenomenon we decomposed 
the calculated binding-energy curves using an empirical formula as suggested in Ref. [23j. In a contrast to standard 
DFT approximations, where asymptotic behavior of the curves cannot be described due to the built in locality [23j, 
in the case of the approach presented in this work, E{d) = Aexjp{—Bd) — C/{d — D)^ represents the perfect fit for 
the computed data-points [see Fig. [s] (Upper panel) as well as Supplementary Fig. SI]. The corresponding energy 
decomposition into repulsive and attractive potential is shown in Fig. [s] (Lower panel). One can see, that indeed, 
Pauli repulsion is weaker, when Xe occupies the fee hollow site. At the same time, an attractive part of the potential 
is weaker in this case as well, when compared to that for the atop geometry. 

As a further step in the understanding, we can explore the advantage of the method of increments to analyze the 
individual contributions to E"^^ computed according to eqn. ([2|. The different contributions to the total adsorption 
energy as obtained for the two considered situation are presented in Fig. |4] As expected, the major contribution 
to the adsorption energy is from the two-body correlation energy increments rjAi- Among them, the interactions 
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between Xe atom and interface Mg layer, (S), are the most important, in particular, one third is due to the electron 
correlations of the adsorbate and the top Mg atom, when the atop geometry is considered [indicated with "TOP" in 
Fig. [4] (Upper panel)]. rjAi involving electrons of Mg-atoms of lower lying layers, (S-1) and (S-2), contribute to the 
^r]Ai with ca. 15% and ca. 2%, respectively. One may note, that the major difference between the two adsorption 
geometries is observed for the increments involving orbitals of entire adsorbate or entire substrate. Their dependance 
from the Xe-Mg(OOOl) distance is plotted in Fig. [s] The curve related to the polarization of the Xe atom is more 
attractive for the atop geometry as compared to the case when adatom is occupying the fee hollow site, indicating 
the higher dipole moment of Xe induced by its interaction with metal surface. Still, the atop- and /cc-related 
curves show similar behavior. This is not the case, when analyzing the results obtained for the sum of the one- and 
two-body increments involving the orbitals of Mg-surface. Whereas the function obtained for the atop geometry is 
attractive for all considered distances, in the case when Xe is occupying the fee hollow site it is weakly repulsive for 
short distances and takes negative values only close to the equilibrium substrate- adsorbate distance. Comparing the 
individual r]ij contributions as obtained for the two considered adsorption geometries, one notes that while in each 
case the correlation contribution is more attractive when Xe is occupying the atop position, the strongest deviation 
can be observed for 7^12, 7714, and 7716 which relate to interaction between orbitals of the (S) and (S-1) Mg-layers (see 
Supplementary Fig. S4). This behavior can be explained in a following way: for the atop geometry, the perturbation 
caused by Xe approaching the Mg surface is effectively screened by the metal (via redistribution of electron density); 
in the case of the fee geometry, such an effect does not take place. Similar conclusion can be done, when comparing 
the difference electron density plots as obtained by means of LDA calculations [23] . 

IV. CONCLUSIONS 

Combining the incremental approach with the special embedding scheme, it is possible to consider metallic systems. 
While DFT with local and semi- local functionals can give reasonable results for cohesive energies, the errors are larger 
for surface stability and adsorption energies, and there is no systematic way of improving the accuracy. On the other 
hand, wavefunction-based correlation techniques, like CCSD(T), allow to get systematically improvable results for 
the full set of characteristics. Thus, the data obtained by means of presented approach can serve as benchmarks in 
the development of improved DFT-based methods. Furthermore, the method of increments allows one not only to 
obtain reference values but also to understand influence of individual correlation-energy increments on the binding 
properties. 

Acknowledgements 



The author would like to thank Prof. Dr. B. Paulus (Berlin) and Prof. Dr. H. Stoll (Stuttgart) for valuable 
discussions. The computing facilities (ZEDAT) of the FU Berlin are acknowledged for computer time. The support 



7 

from the German Research Foundation (DFG) through the Cohaborative Research Center (SFB) 765 is appreciated. 
Worthful suggestions and useful discussions during preparation of the manuscript from Dr. Yu. S. Dedkov (Berhn) 
are acknowledged. 



[1] P. Pulay, Chem. Phys. Lett. 100, 151 (1983). 

[2] S. Saeb0 and P. Pulay, Chem. Phys. Lett. 113, 13 (1985); P. Pulay and S. Saeb0, Theor. Chim. Acta 69, 357 (1986); S. 

Saeb0 and P. Pulay, J. Chem. Phys. 86, 914 (1987); S. Saeb0 and P. Pulay, J. Chem. Phys. 88, 1884 (1988). 
[3] M. Schiitz, G. Hetzer, and H.-J. Werner, J. Chem. Phys. Ill, 5691 (1999). 
[4] C. Hampel and H.-J. Werner, J. Chem. Phys. 104, 6286 (1996). 
[5] M. Schiitz and H.-J. Werner, J. Chem. Phys. 114, 661 (2001). 

[6] M. Schiitz and H.-J. Werner, Chem. Phys. Lett. 318, 370 (2000); M. Schiitz, J. Chem. Phys. 113, 9986 (2000); M. Schiitz, 

J. Chem. Phys. 116, 8772 (2002). 
[7] H. Stoh, Phys. Rev. B 46, 6700 (1992). 
[8] B. Paulus, Phys. Rep. 428, 1 (2006). 

[9] E. Voloshina, N. Gaston, and B. Paulus, J. Chem. Phys. 126, 134115 (2007). 
[10] B. Paulus and K. Rosciszewski, Chem. Phys. Lett. 394, 96 (2004); B. Paulus, K. Rosciszewski, N. Gaston, P. Schwerdtfeger, 

and H. StoU, Phys. Rev. B 70, 165106 (2004); N. Gaston, B. Paulus, K. Rosciszewski, P. Schwerdtfeger, and H. Stoh, 

Phys. Rev. B 74, 094102 (2006). 
[11] E. Voloshina and B. Paulus, Phys. Rev. B 75, 245117 (2007); Mol. Phys. 19-22, 2849 (2007). 

[12] N. Gaston and B. Paulus, Phys. Rev. B 76, 214116 (2007); N. Gaston, B. Paulus, U. Wedig, and M. Jansen, Phys. Rev. 
Lett. 100, 226404 (2008). 

[13] N. Gaston, D. Andrae, B. Paulus, U. Wedig, and M. Jansen, Phys. Chem. Chem. Phys. 12, 681 (2010). 

[14] E. Voloshina and B. Paulus, Z. Phys. Chem. 224, 369 (2010); E. Voloshina and B. Paulus, Chem. ModeU. 6, 162 (2009). 

[15] B. Paulus, Surf. Sci. 408, 195 (1998). 

[16] E. Voloshina, D. Usvyat, M. Schiitz, Yu. S. Dedkov, and B. Paulus, Phys. Chem. Chem. Phys. 13, 12041 (2011). 

[17] E. Voloshina and B. Paulus, Theor. Chem. Acc. (2005) 114,259 (2005); E. Voloshina and B. Paulus, J. Chem. Phys. 124, 

234711 (2006); E. Voloshina and B. Paulus, J. Comput. Chem. 29, 2107 (2008). 
[18] C. Miiller, B. Herschend, K. Hermansson, and B. Paulus, J. Chem. Phys., 128, 214701 (2008); C. Miiller, B. Paulus, and 

K. Hermansson, Surf. Sci. 603, 2619 (2009); C.Miiller, K. Hermansson, and B.Paulus, Chem. Phys. 362, 91 (2009). 
[19] B. Paulus and K. Rosciszewski, Int. J. Quantum Chem. 109, 3055 (2009). 
[20] W. R. Tyson and W. A. Miller, Surf. Sci. 62, 267 (1977). 

[21] Th. Seyller, M. Caragiu, R.D. Diehl, P. Kaukasoina, M. Lindroos, Chem. Phys. Lett. 291, 567 (1998); Th. Seyller, M. 

Caragiu, R. D. Diehl, P. Kaukasoina, and M. Lindroos, Phys. Rev. B 60, 11084 (1999). 
[22] P. Thomas, J. Gray, X.D. Zhu, and C.Y. Pong, Chem. Phys. Lett. 381, 376 (2003). 

[23] J. L. F. Da Silva, C. Stampfl, and M. Scheffler, Phys, Rev. Lett. 90, 066104 (2003); J. L. F. Da Silva, C. Stampfl, and M. 



Scheffler, Phys, Rev. B. 72, 075424 (2005); J. L. F. Da Silva, Ph.D. thesis ( |http://www.fhi-berlin.mpg.de/th/[ pub02.html). 



8 



Technical University Berlin, Berlin, Germany, 2002. 
[24] P. A. M. Dirac, Proc. Camb. Phil. Soc. 26, 376 (1930); S.H. Vosko, L. Wilk, and M. Nusair, Can. J. Phys. 58, 1200 (1980). 
[25] J. P. Perdew, K. Burke, and M. Ernzerhof, Phys. Rev. Lett. 77, 3865 (1996). 

[26] A. Stroppa and G. Kresse, New. J. Phys. 10, 063020 (2008); L. Schimka, J. Harl, A. Stroppa, A. Griineis, M. Marsman, 

F. Mittendorfer, and G. Kresse, Nature Materials 9, 741 (2010). 
[27] S. Grimme, J. Comput. Chem. 27, 1787 (2006). 

[28] R. Dovesi, R. Orlando, B. Civalleri, C. Roetti, V.R. Saunders, CM. Zicovich- Wilson, Z. Kristallogr. 220, 571 (2005); R. 

Dovesi, V. R. Saunders, C. Roetti et al, CRYSTAL09 Users Manual, University of Torino, Torino, 2009. 
[29] U. Wedig, M. Jansen, B. Paulus, K. Rosciszewski, and P. Sony, Phys. Rev. B 75, 205123 (2007). 
[30] A. Nicklass, M. Dolg, H. Stoll, H. Preuss, J. Chem. Phys. 102, 8942 (1995). 

[31] MOLPRO, version 2006.1, a package of ab initio programs, H.-J. Werner, P. J. Knowles, F. R. Manby, M. Schiitz, and 



others, see http://www.molpro.net 



[32] C. Hampel, K. Peterson, H.-J. Werner, Chem. Phys. Lett. 190, 1 (1992); J. D. Watts, J. Gauss, R. J. Bartlett, J. Chem. 

Phys. 98, 8718 (1995); M. J. O. Deegan, P. J. Knowles, Chem. Phys. Lett. 227, 321 (1994). 
[33] T.H. Dunning, Jr. (unpublished). 

[34] J. M. L. Martin and A. Sundermann, J. Chem. Phys. 114, 3408 (2001). 
[35] P. Fuentealba et a/., J. Phys. B 18, 1287 (1985). 

[36] C. Kittel, Introduction to Solid State Physics, 7th ed., (John Wiley and Sons, New York, Chichester, 1996). 



9 



TABLE I: Lattice constants (a, c), cohesive ), surface energy (^surf), equilibrium distance (do), and adsorption 

energy (^ads) evaluated for the hep Mg, Mg(OOOl) surface, and Xe/Mg(0001) as computed for two adsorption geometries (c./. 
text), respectively, as obtained at different levels of theory. 



Method 




hep Mg 


Mg(OOOl) 


Xe/Mg(0001) 


atop position 


fee hollow site 


a (A) 


c(A) 


^coh (eV/atom) 


^surf (eV/atom) 


do (A) 


^ads (meV) 


do (A) 


^ads (meV) 


Expt. 


3.21"^ 


5.21"^ 


-1.51" 


+0.383^ 










HF 


3.31" 


5.13" 


-0.36" 


+0.182 


no binding 




no binding 




CCSD(T)^ 


3.21" 


5.14" 


-1.52" 


+0.385 


4.12 


-72 


4.28 


-67 


LDA 


3.13" 


5.00" 


-1.80" 


+0.404 


3.83 


-91 


3.87 


-85 


PBE 


3.19" 


5.12" 


-1.48" 


+0.360 


4.81 


-13 


4.94 


-12 


PBE-D2 










4.04 


-75 


4.13 


-70 



"Ref. 1361: ^Ref. |20 | : "^Ref. [29 | : ^Correlation contribution to E^r^u, E^.^yf, and E^r^^ are obtained according to eans. ([l[). 1^, and ([2] 
respectively; ^Ref. 



10 



TABLE II: Individual 2-body correlation energy increments (in eV) as obtained for "bulk" and "surface" clusters with CCSD(T) 
method. Correlation contributions to the surface energy (7) are obtained when using eqn.js] (Note: 1- and 3-body contributions 
are listed in Supplementary Tab. SI. Structure-related data for the considered increments are summarized in Supplementary 
Figs. S2 and S3). 









, hulk 




surf 




12 


3.197 


-0.15411 


3 


-0.18050 


3/2 


+0.19158 


13" 


3.209 


-0.15190 


3 


-0.17156 


3 


-0.05898 


14 


4.529 


-0.02768 


3 


-0.03299 


3/2 


+0.03355 


15 


5.210 


-0.01375 


1 


-0.02217 


1/2 


+0.00267 


16 


5.551 


-0.00982 


6 


-0.01272 


3 


+0.02079 


17" 


5.558 


-0.00970 


3 


-0.01320 


3 


-0.01051 


18 


6.119 


-0.00591 


6 


-0.00971 


3 


+0.00633 


19" 


6.418 


-0.00515 


3 


-0.00730 


3 


-0.00645 



"the both atoms of the central part of the cluster belong to the same layer of the hep lattice 
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o 



(S) 



(S-1) 



(S-2) 




Central part: 
Q VTZ 



Embedding: 

• Mg2^-PP + VDZ 

e (o) Mg2+-PP + minimal basis 



FIG. 1: Example of a "bulk" Mg-cluster for calculating the 2-body contribution to the correlation energy (ri2 = 3.197 A). 
Spheres of different size and style represent atoms of different types (c./. text). For the one of the central atoms, six nearest- 
neighbor distances are depicted with thick lines (either solid or dashed). Atoms shown with plain spheres (either hatched or 
filled) disappear when constructing the "surface" cluster. Crystallographic structure of hep Mg is shown with thin dashed lines. 
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FIG. 2: Correlation contributions to the total energy of 1-, 2-, and 3-body increments and calculated surface energy as 
compared with experimental value. 
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FIG. 3: (Upper panel) The different contributions to the adsorption energy as functions of substrate- adsorbate distance are 
plotted for the Xe/Mg(0001) system as obtained for the atop adsorption geometry. Circles represent the data points, obtained 
by summing up the HF (thick solid line) and correlation (dashed line) parts of energy. Thin solid line is a fit according to the 
formula E(d) = Ae:s.p(—Bd) — C/{d — D)^ [c.f. text). (Lower panel) Decomposition of the binding energy curve into repulsive 
and attractive contributions {c.f. text). 





FIG. 4: The different contributions to the total adsorption energy are plotted for the Xe/Mg(0001) system in its equilibrium 
geometry, when Xe is occupying the atop position (upper panel) or the fee hollow site (lower panel). 
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FIG. 5: Correlation contributions to the adsorption within the adsorbate (?7a) and within the surface C^r/i -\- ^r/ij) as 
functions of distance between Xe and the magnesium surface. Equihbrium distances obtained for the two adsorption positions 
are shown with vertical lines. 
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Supplementary material for the manuscript: 

Local correlation method for metals: benchmarks for surface and adsorption energies 

Elena Voloshina 

Physikalische und Theoretische Chemie, Freie Universitdt Berlin, 14195 Berlin, Germany 
Content: 

1. Tab. SI: Individual correlation- energy increments (in eV) as obtained for "bulk" and "surface" clusters with 
CCSD(T) method. 

2. Fig. SI: The different contributions to the adsorption energy as functions of substrate- adsorbate distance are 
plotted for the Xe/Mg(0001) system as obtained for the two different adsorption geometries. Circles represent 
the data points, obtained by summing up the HF (thick solid line) and correlation (dashed line) parts of energy. 
Thin solid line is a fit according to the formula E{d) = Aexp{—Bd) — C/{d — D)^ {c.f. text). 

3. Fig. S2: The two-body correlation-energy increments are plotted as a function of distance. 

4. Fig. S3: The considered 3-body structures are shown schematically on the background of a hexagon in the hep 
plane. 

5. Fig. S4: Correlation contributions to the adsorption energy within the substrate as functions of the substrate- 
adsorbate distance, when considering the Xe/Mg(0001) system as obtained for (a) one-body increment r]i; (b-d) 
two body increments r]ij. (b): 7/13, 7717, 7/19 indicate interaction between orbitals within the (S) layer; (c): 7712, 
^14, V16 indicate interaction between orbitals of the (S) and (S-1) layer; (d): 7/15, rjig indicate interaction between 
orbitals of the (S) and (S-2) layer. The energies are weighted for the total contribution of each increment. The 
scale is not the same in each subplot. 
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Table SI. Individual correlation-energy increments (in eV) as obtained for "bulk" and "surface" clusters with CCSD(T) 
method. Correlation contributions to the surface energy (7) are determi- ned as a difference between the increment calculated 
using the "surface" cluster, multiplied with its weight factor (uj^^^^), and the corresponding one calculated using the "bulk" 
cluster, multiplied with its weight factor (cj^^^^). 





"bulk" 


^bulk 


"surface" 


^surf 


7 


1 


-0.59990 


1 


-0.65230 


1 


-0.05240 


E l-body 


-0.59990 




-0.65230 




-0.05240 


12 


-0.15411 


3 


-0.18050 


3/2 


+0.19158 


13" 


-0.15190 


3 


-0.17156 


3 


-0.05898 


14 


-0.02768 


3 


-0.03299 


3/2 


+0.03355 


15 


-0.01375 


1 


-0.02217 


1/2 


+0.00267 


16 


-0.00982 


6 


-0.01272 


3 


+0.02079 


17" 


-0.00970 


3 


-0.01320 


3 


-0.01051 


18 


-0.00591 


6 


-0.00971 


3 


+0.00633 


19" 


-0.00515 


3 


-0.00730 


3 


-0.00645 


E 2-body 


-1.15375 




-0.97479 




+0.17896 


3.1[2+1]^ 


-0.00276 


3 


+0.00785 


5/2 


+0.02791 


3.1[l+2]^ 


-0.00276 


3 


+0.01637 


2 


+0.04102 


3.2" 


+0.00198 


1 


+0.00534 


1 


+0.00336 


3.3" 


-0.00596 


1 


-0.00272 


1 


+0.00324 


3.4[2+l]^ 


+0.00078 


6 


+0.00173 


5 


+0.00397 


3.4[l+2]'' 


+0.00078 


6 


+0.00756 


4 


+0.02556 


3.5 


-0.00134 


3 


-0.00324 


5/2 


-0.00408 


3.6[2+l]^ 


-0.00132 


6 


-0.00498 


5 


-0.01698 


3.6[l+2]^ 


-0.00132 


6 


+0.00076 


4 


+0.01096 


3.7" 


-0.00130 


6 


-0.00163 


6 


-0.00198 


3.8 


-0.00244 


6 


-0.00582 


5 


-0.01446 


3.9" 


-0.00271 


3 


-0.00320 


3 


-0.00147 


E 3-body 


-0.06161 




+0.01544 




+0.07705 


TTicorr 
-t^total 


-1.81526 




-1.61165 




+0.20361 



"all atoms of the central part of the cluster belong to the same layer (S) of the hep lattice; 

^two atoms of the central part belong to the (S)-plane and the remaining atom is lying in the (S-1) plane; 

*^one atom of the central part belongs to the (S)-plane and the remaining two atoms are lying in the (S-1) plane 
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Fig. SI. The different contributions to the adsorption energy as functions of substrate- adsorbate distance are plotted 
for the Xe/Mg(0001) system as obtained for the two different adsorption geometries. Circles represent the data points, 
obtained by summing up the HF (thick solid line) and correlation (dashed line) parts of energy. Thin solid line is a 
fit according to the formula E{d) = Aeyip{—Bd) — C/{d — D)^ {c.f. text). 




Fig. S2. The two-body correlation-energy increments are plotted as a function of distance. Inset: In the crystallo- 
graphic structure of hep Mg, the atoms are ordered by their distance from the atom number 1. 
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Fig. S3. The considered 3-body structures are shown schematicahy on the background of a hexagon in the hep plane. 
Light spheres correspond to atoms below the plane, and crosses to atoms above the plane. 
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Fig. S4. Correlation contributions to the adsorption energy within the substrate as functions of the substrate- 
adsorbate distance, when considering the Xe/Mg(0001) system as obtained for (a) one-body increment 77^; (b-d) two 
body increments r]ij. (b): 7713, 7717, 7^19 indicate interaction between orbitals within the (S) layer; (c): 7712, 7714, 7^16 
indicate interaction between orbitals of the (S) and (S-1) layer; (d): 7^15, r]is indicate interaction between orbitals of 
the (S) and (S-2) layer. The energies are weighted for the total contribution of each increment. The scale is not the 
same in each subplot. 



